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Q> ' The running of the couphng is studied in SU(4) gauge theory using the Schrodinger 

■ functional technique. Up to energies of the order of the square root of the string tension 

I— —I, a, the running is found to agree with the two-loop perturbative formula. Relating 

I the perturbative to the non-perturbative regime of the running and converting to the 

I ■ MS scheme allows one to extract the ratio Aj^/ ^/a. The result is then used in 

combination with similar calculations present in the literature for SU(2) and SU(3) 
r~| \ to extract Ajj^/^/a in the large limit. Our results for N = 3,4 agree with a 

recent study of the same quantity performed using the Parisi mean field scheme as an 
^ . intermediate scheme, while K-^/y/a in SU(2) turns out to differ by 2.5%. Possible 

^ \ explanations of this discrepancy are discussed. 

(N ; 1 Introduction 

in ■ 

I Asymptotic freedom [1,2] is one of the signatures of non-Abelian gauge theories. At 

short distances, the coupling is small and perturbation theory can be successfully used 
to compute observables. By dimensional transmutation, a scale is generated (of mass 
.5^ \ dimension 1) whose order of magnitude is dictated by the dynamics of the theory. 

^ ' In perturbation theory, it is convenient to associate with this scale the dimensionful 

^ . multiplicative constant of the integrated perturbative beta function, which is called 

the A parameter. The value of the A parameter depends on the chosen perturbative 
scheme. Different regularisation schemes define different couplings. A good scheme is 
conventionally one for which low-order perturbative calculations work at energies close 
to the non-perturbative scale of the theory. 

However good the scheme is, at long enough distances (low energies) the theory 
becomes inherently non-perturbative, and confinement sets in. In SU(iV) gauge theo- 
ries (or more generally in theories with unbroken fundamental centre), this regime is 
characterised by the tension of the confining string, cr, which has dimensions of mass 
squared. Since SU(A^) gauge theories are theories with only one dynamically generated 
scale, it must be possible to relate A and a. For SU(2) and SU(3) gauge theory, this 
program has been successfully carried out by the Alpha collaboration [3-5] using the 
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Schrodinger functional (SF) scheme [6, 7] on a spacetime lattice. In the SF scheme, 
the theory is defined in a box of finite physical size, and the renormalised coupling 
is obtained through the effective action of the system with certain specified boundary 
conditions. This determines the coupling as a function of the extension of the box, 
L. On the lattice, simulations are carried out at different lattice spacings in order to 
extract continuum results. In doing so, the bare coupling must be carefully tuned for 
each lattice spacing to ensure that the physical size L remains constant. An iterative 
procedure can then be set up to probe the theory over a large range of energies. 

This technique has proved to be effective for exploring a range of couplings interpo- 
lating from the perturbative to the non-perturbative regime. On the perturbative end, 
the A parameter in the SF scheme, A^f, can be determined using perturbation theory, 
in terms of the physical size of the box; on the non-perturbative end, the size of the box 
can be determined in terms of a non-perturbative quantity, e.g. the string tension a. 
Since the relationship between the sizes of the box at the two ends is known, one can 
determine the ratio Asf / v^- By calculating Asf / using perturbation theory, the 
ratio can be re-expressed in the language of the more familiar MS scheme as Ajj^/^/a. 

In recent years, following renewed interest from string theory [8], a program for 
non-perturbative studies of SU(A^) gauge theories in the large-A'' limit [9] has been 
developed (see e.g. [10-21]). A general conclusion is that observables in SU(A^) pure 
gauge theories have a smooth dependence on N that, within a few percent, can be 
accounted for by a 1/iV^ correction to the N = oo case for A'^ > 3, and often also 
including the case N = 2. In this work, we investigate the dependence of Ajj^/^/a as 
a function of A^ using the SF technique. To this end, first we formulate the problem 
of the running of the coupling for SU(N) gauge theory in the SF scheme (Sect. [2]) 
and numerically determine the value of Aj^/ ^/a for SU(4) following the procedure 
described in [3,4] (Sect. [3l). Then, using results available in the literature for SU(2) 
and SU(3) complemented with recent lattice determinations of the string tension, we 
discuss the behaviour of Ajjg/^/a as N ^ oo (Sect. [4]), comparing our results with 
those of [22], where the Parisi mean field improvement [23] is used as an intermediate 
scheme. A discussion of systematic errors follows (Sect. [5]). Finally, we draw our 
conclusions in Sect. El A partial account of our calculation has already been published 



2 The Schrodinger functional in SU(A^) gauge 
theory 



Following Ref. [7] (to which we refer for further details), we shall introduce the SF in 
SU(A'") lattice gauge theory. Consider two states C and C in the Schrodinger represen- 
tation of a system described by the Hamiltonian H, whose associated action is given 
by 5. The Schrodinger functional is the probability amplitude for C at time t starting 
from C at time 0: 



in [24]. 
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where 4> is the generic field configuration and D[(f)] the measure of the path integral, 
which is taken at fixed boundary conditions. For a SU(A'^) lattice gauge theory on a 
volume V = with L = na {n being an integer and a the lattice spacing), described 
by the Wilson action, Eq. ([T]) becomes 



Z[C,C'] 



D[U]e 



-S[U] 



-r[c,c'] 



(2) 



C,C' 



where r[C,C] is the effective action of the system with the specified boundary condi- 
tions. The Wilson action appearing in the above definition is given by 



S[[/] = l^Tr(l-t/(p)), 



(3) 



where U{p) denotes the parallel transport of the link variables U^{x) (x = {x^ , x^, x^, x^) 
being a lattice point of integer coordinates and /i = 0, 1, 2, 3 the lattice directions) over 
the elementary square of the lattice (plaquette) p and go is the (bare) lattice coupling. 
The sum in Eq. ([3]) must be taken over both orientations of the plaquettes. 

The boundary links W are required to satisfy inhomogeneous Dirichlet boundary 
conditions for /c = 1, 2, 3, 



Wk{x)\^o=o = exp (aCfc(x)) , Wk{x)\^o=L = exp (aCfc(x)) 



(4) 



where Ck and are spatial boundary fields which need to be chosen. 

Fixing the boundary field induces a background field in the bulk, and it is desirable 
to choose boundary fields which minimise the effect of the finite lattice spacing. It was 
shown in [7] that, for A'^ colours, the optimal choice are constant Abelian fields. 
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Unitarity and stability considerations of the background field constrain the angles, 

N 



i=l 



0, 



^l<4'2 < ■■■ < (pN, 



< 2tt 



(6) 



and similarly for (from here on we drop the suffix k on the angles and use the same 
choice for k = 1,2,3). The ensemble of points satisfying the constraints ^ is referred 
to as the fundamental domain. 

The effective action can be written as an asymptotic series 



T[B] = go'TolB] + Ti[B] + giV^lB] + 



(7) 



In this expansion, Fq is the classical action, which can be evaluated analytically. The 
previous equation could be used directly to define a renormalised coupling via Monte 
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Carlo simulations of the SF. However, the numerical determination of an effective action 
is a notoriously difficult problem. To get around this difficulty, one generally measures 
derivatives of the effective action. By introducing a dependence of the boundary links 
(and thus of the background field) on a real dimensionless parameter rj, we can then 
define a renormalised coupling as 

f = ^^. r'iBi = |-r[Bi. (8) 

for a particular choice of rj. 

The SF technique allows one to explore the running of the coupling for a wide 
range of energies, connecting the perturbative to the non-perturbative regime. On the 
lattice, for each energy scale Ei probed, in order to extrapolate to the continuum limit, 
multiple simulations must be performed at different lattice spacings; hence, different 
bare couplings are needed. Those couplings (ordered from the one corresponding to the 
coarsest to the one corresponding to the finest simulated lattice spacing) are labeled 
sequentially by a second index j. Thus, we define g^-' as the jth bare coupling at 
the energy scale Ei, or, equivalently, the length scale Lj = 1/Ei. The running of the 
coupling is computed using a recursive procedure. We start by fixing the renormalised 
coupling for all ^q'-^ in such a way that it is equal for the whole set, within errors. 
This common value of the renormalised coupling is called g'^. Then, for each bare 
coupling 5(q'"' , the renormalised coupling is evaluated for boxes of size 2Lq (by doubling 
the number of lattice sites in each direction) and the continuum limit is obtained 
by assuming a linear dependency in a/L and extrapolating to a ^ 0. We call the 
extrapolated value g^. Returning to a small lattice once again, a new set of couplings 
gQ'' is then chosen in such a way that the value of the renormalised couplings for the 
size Li match g^. This ensures that the physical size of the box for the set g^'-^ is kept 
constant, and in particular equal to twice the size of the box corresponding to the set 
^q'-'. This procedure can then be iterated. If Lq is the original size of the box, after I 
iterations the size is Li = 2^Lq. If the final set of couplings g^ are in the asymptotic 
scaling regime of the theory, the product Li\fo can be obtained by determining a for 



each coupling in the set g^K On the other end, if g^ is in the perturbative regime, Lq 
can be obtained from the integrated two-loop beta function of the theory as 



where 



Putting together the determination of = Lmax in terms of ^Jo and h^gp^ the value 
of KsF/\f^ can be worked out. 

In order to compare with other determinations, the SF scheme should be related to 
a more widely used scheme, e.g. the MS. Two regularisation schemes can be related 
by using first order perturbation theory. Consider two schemes, A and B, in which the 
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couplings QA and qb are defined. The corresponding A parameters are found through 
the one-loop beta functions: 

1 , 1 



K = — e ^i^o'A and K = —e ^''o^s , (11) 
Aa As 

where K is the length scale at which the coupling is measured. Putting together the 
two previous relationships yields 

^=e"^fe"4). (12) 
Ab 

Prom a first order perturbative calculation of a physical quantity one gets 

9a 9% 

where A; is a constant depending on the details of the schemes. Once k has been 
determined, we can rewrite the ratio of the A parameters as 

^ = e~^^ . (14) 
As 



To relate the SF and MS scheme we can use the lattice scheme as an intermediate 
step, computing Asf/^l and Aj^^^/A^. The perturbative calculation relating the 
lattice and the SF schemes is described in [7], while the ratio Aj^/Aj^ for SU(A'') 
gauge theories (determined in Ref. [25]) is given by 

= 38.85e"n^ . (15) 

Al 

Other perturbative calculations of A-j^/A^ [26,27] give slightly discrepant results. 
However, the uncertainty coming from this perturbative calculation is some orders of 
magnitude smaller than the error on Asf/^l, and can safely be neglected. 



3 The running of the couphng in SU(4) gauge 
theory 

We now specialise to SU(4). The fundamental domain can be described symmetrically 
by defining a one-to-one map between the set of angles (^i, <?!>2, ^/'s, ^4) and a point V 
in a certain bounded three-dimensional region, 

V = 1(^1 • ei -1-^2 • 62 -|-<;/;)3 ■ 63 -1-^4 • 64), <pi = Y-ei, (16) 

where are the weights of the Lie algebra of SU(4) in the fundamental representation, 
normalised as • = |(4(5ij — 1). With this normalisation, the vertices of the 
fundamental domain are 

Vi = -Y^e,-, z = 1,2,3,4, (17) 
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Figure 1: The fundamental domain and its two planes of symmetry. The points marked 
correspond to the angles (fTSi) . 

describing a skewed tetrahedron (Figure [1]). One must then select, from inside the 
fundamental domain, the set of angles for each boundary, and C^, each set cor- 
responding to a point inside the tetrahedron. These points will be members of a 
one-parameter family of angles parameterised by rj, through which the renormalised 
coupling is defined by ([8|). In principle this choice is arbitrary and of no conceptual 
significance; however the signal-to-noise ratio of the Monte Carlo evaluation is highest 
when (i) the points are well away from the domain's edge, (ii) they are as far as possible 
from each other, and (iii) the two boundaries are on an equal footing. Geometrically, 
this corresponds to choosing two points related by a symmetry of the fundamental 
domain. 

We choose the two points to be related by the symmetry refiecting about the plane 
through V2 and V4, and make the particular choice 



h = -hv-hV2 </>'i = ir/-i(2 + ^/2)7r 

^2 = -h - iC^ - (/.'2 = ir/-i(4-^/2)7r 

(/>3 = i7? + i(2-^/2)7r 0?j = -ir? + i(4- V2)7r 

(j)4 = h + h^ </.^ = -i77 + i(2 + ^/2)7r 



With this choice, we set rj = and compute the renormalised coupling ^ by calculating 



the expectation value of the observable 



1^ = E E [(^^ w + w) + w + w)'' 




where c = diag(— ^, — ^, ^, ^), and a similar expression holds for i?'(x). In our simu- 
lations, we have used the Cabibbo-Marinari algorithm, with one heat-bath update for 
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P = 


L/a -g\L) g\2L) 




15.126 
15.626 
16.000 


6 1.0222(6) 1.2247(10) 
8 1.0223(7) 1.2162(18) 
10 1.0223(5) 1.2104(21) 


1.1892(11) 


14.137 
14.632 
15.007 


6 1.1893(3) 1.4833(13) 
8 1.1892(4) 1.4705(19) 
10 1.1890(8) 1.4635(25) 


1.4331(18) 


13.142 
13.631 
14.000 


6 1.4329(3) 1.9043(14) 
8 1.4332(4) 1.8798(24) 
10 1.4331(6) 1.8668(35) 


1.8099(24) 


12.190 
12.668 
13.030 


6 1.8098(7) 2.7198(38) 
8 1.8102(6) 2.6548(49) 
10 1.8094(6) 2.6177(49) 


2.4645(33) 



Table 1: Pairs of renormalised couplings for fixed values of f3 = S/qq. 

every four over-relaxation steps. For the highest values of (3 = S/g^ (i.e. smallest bare 
coupling), the number of configurations generated for each data point is of order 10^, 
and this number increases for decreasing (3, up to ~ 10^ configurations for the lowest 
p. 

Following the outline of Sect. [21 we have performed four sets of simulations, each 
at three different lattice spacings. The results are shown in Tab. [TJ Using the data 
in columns 2 and 4, we have extrapolated the renormalised coupling to the continuum 
using a linear fit in a/L, with a reduced of ~ 0.1 for each fit. Although in principle 
higher order corrections in a/L would need to be considered, the quality of the fit 
(showed by the low value of the x^) suggests that the dominating contribution comes 
from the term proportional to a/L. The extrapolation is plotted in Fig. [2] and the result 
is shown in column 5 of Tab.[TJ All the errors shown are statistical. Note, however, that 
there is an accumulation of errors due to the iterative nature of the procedure. This 
is because the statistical error in the extrapolated value of the renormalised coupling 
for some set i (column 5) leads to an uncertainty on the value to which we tune the 
renormalised coupling in the next set z + 1 (column 3). As in [3-5], it is convenient 
to interpret this mismatch as an error on the scale going from Lj to Lj+i, and thus 
also from Lq to Lmax- The mismatch is small, so we can integrate the two-loop beta 
function to estimate the implied error on the scale at each stage, and then add the 
errors together in quadrature (since they come from independent sets of simulations). 
We find that Lmax = (2^ ± 0.04)Lo. 

At the most non-perturbative point (i.e. for = 2.4645), we can express the 
scale in units of the string tension. In Tab. [21 we quote the bare coupling at fixed 
renormalised coupling for different values of L/a, and the corresponding string tension 
^/a in lattice units extracted using the interpolating formulae in [14] supplemented by 
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Figure 2: Extrapolation to the continuum of the couphngs in Tab. [TJ 



L/a 


/9 


9\L) 




LhiAX^f^ 


5 


11.029 


2.4644(4) 


0.2093(10) 


1.046(5) 


6 


11.326 


2.4646(7) 


0.1611(10) 


0.966(6) 


7 


11.574 


2.4645(18) 


0.1340(10) 


0.938(7) 


8 


11.782 


2.4646(8) 


0.1146(10) 


0.917(8) 



Table 2: Data used for the extrapolation of LMAx\f^ to zero lattice spacing for SU(4) gauge 
theory. 
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Lja 






LMAX\f^ 


5 


2.5009 


0.1835(10) 


0.918(5) 


6 


2.5752 


0.1438(10) 


0.863(6) 


7 


2.6376 


0.1179(10) 


0.825(7) 


8 


2.6957 


0.0984(10) 


0.787(8) 


9 


2.7378 


0.0859(10) 


0.773(9) 


10 


2.7824 


0.0751(10) 


0.751(10) 



Table 3: Data used for the extrapolation of Lmax\/^ to zero lattice spacing for SU(2) gauge 
theory. The j3 values have been taken from [3] and the string tensions from [14,22]. 

newer data in [220. The values ol Lmax thus obtained then need to be extrapolated to 
the continuum. This requires a delicate fit, as it has been found (see e.g. the discussion 
in [28]) that lattice corrections are large in both leading and next-to- leading order. We 
fit using 

LmAX\/o- = LMAxV^\a=0 + C\a^ + 020^0 , (20) 

and estimate the systematic error by varying the extremes of the fit and noting how 
much the extrapolated value shifts. This leads to a value oiLMAX\f^ = 0.910(40)(100) 
where the first parentheses give the statistical error of the fit, and the second give our 
estimate of the systematic error. Using Q, this gives 

AsF = 0.253(10) (30) Vo^= 106.3 MeV ±4.2 MeV ±12.6 MeV . (21) 

where the value in MeV is obtained using yfa = 420 MeV. Finally, we have calculated 
the ratio 

hrMsl^SF = 2.08114(34) (22) 

and, at the same time, calculated the improvement coefficient c^^^ = —0.12005(15) (in 
the notation of Ref. [7] ) . This immediately leads to 

^M5/^^ = 0-527(21)(62) . (23) 

We plot the renormalised coupling against the energy scale using the lattice data, 
together with the one- and two-loop perturbative predictions in Fig. [3l As was found 
in the previous studies of SU(2) and SU(3), two- loop perturbation theory gives excellent 
agreement with the data up to the scale of the string tension. 

4 Large- limit 

Using the data for the SF from Refs. [3, 4] and the data for the string tension from 
Refs. [14,22], we can compute the ratio K-JJg/^/a for SU(2) and SU(3) gauge theory 

^The results a,i (3 — 11.574 and /? = 11.782 were obtained by extrapolating an unpublished string tension 
measurement at (3 — 11.5. 
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Figure 3: Results of lattice simulations for the running coupling in the SU(4) theory, together 
with the one- and two-loop perturbative predictions. Energies in physical units are obtained 
by noting ^ = 420 MeV. 



L/a 


P 


a^/o' 




4 


5.9044 


0.2559(10) 


1.024(4) 


5 


6.0829 


0.1884(10) 


0.942(5) 


6 


6.2204 


0.1528(10) 


0.917(6) 


7 


6.3443 


0.1282(10) 


0.897(7) 


8 


6.4527 


0.1107(10) 


0.885(8) 



Table 4: Data used for the extrapolation of LuAx^f^ to zero lattice spacing for SU(3) gauge 
theory. The /? values have been taken from [4] and the string tensions from [14,22]. 
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^Ms/ ^SF 


2 


2.2446 


3 


2.0487 


4 


2.0811 



Table 5: Conversion coefficients from the SF scheme to the MS scheme. 

in analogy with the SU(4) case discussed in the previous section. The only difference 
appears in the case of SU(3), where the simulations have been done with a Symanzik- 
improved action. Here we expect the linear coefficient in (I20p to be small (but not zero 
as the improvement is only calculated to 1-loop in perturbation theory), and indeed 
we observe this in our fit where we find ci = —0.20 compared to C2 = 3.12. We also fit 
with ci set to zero, and take the difference as an indication of the systematic error. For 
convenience, we report the data we have used for the extrapolations in Tabs. [3] and [H 
In the continuum we find 

^MAX\/^ = 0.603(17) (50) (24) 

for SU(2) and 

LmaxV^ = 0-8M{3){m) (25) 
for SU(3). Using the conversion factors in Tab. [U those results imply 

AjTs/Vo^ = 0.752(20)(60) (26) 

and 

^Ms/^^ = 0-538(l)(20) (27) 

for SU(2) and SU(3) respectively, where the first parenthesis is the statistical error and 
the second is an estimate of the systematic error coming from the uncertain form of 
the fit. 

A reliable determination of the large- limit of Aj^/^/a with the SF technique 
requires the determination of this quantity at larger values of N, which is outside the 
scope of this work. However, since a calculation in the large- A^ limit of this quantity 
has been provided in [22], we can compare our data with the extrapolation reported 
there, to check the effects of possible systematic errors. We find that the results for 
N = 3,4: are in good agreement, while for SU(2) the value of Ajjg/y/a obtained 
with the SF technique is higher (2.5% when systematic errors of both calculations 
are included). However, there are indications that the SU(2) case is problematic in 
both calculations: with the SF technique, the extrapolation to the continuum limit 
proves to be less controlled than for SU(3) and SU(4) (an issue that can be resolved 
in principle by including larger volumes or using an improved action to determine the 
running of the coupling), while for the Parisi mean field method the plaquette in SU(2) 
might be affected by the end point of the bulk transition in the fundamental-adjoint 
plane (we refer to [22] for further details). Although perfectly reasonable at this level, 
this discrepancy deserves further investigation. The agreement of the = 3, 4 values 
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1/N 



Figure 4: The ratio A-^/ y/a determined in this work as a function of l/N"^. The continuous 
hne is the extrapolation to = cxd from Ref. [22], while the dot-dashed lines delimit the 
region at one sigma of confidence level (only the statistical errors are shown). 

computed with two independent methods shows that at large the systematic should 
be under control in both cases. 



5 Systematic errors 

As we have seen in the previous section, albeit with larger errors, our results agree 
with those reported in [22]. In order to better assess the scope of our findings, in this 
section we shall discuss in detail the systematic errors of our calculation. 

The main sources of systematic errors are the following: 

1. Extrapolation of Lmax\^ to a = 0. We have already mentioned that this is the 
biggest source of error in determining the A-parameter. In order to extrapolate 
LmaxV^ to a = 0, the string tension is measured for a number of lattice spacings 
and then the a — > limit is taken according to Eq. (j20p . The difficulty is the 
largeness of both coefficients ci and C2. We also have no information about the 
rest of the series. By changing the extremes of the fitting interval and compar- 
ing with fits including also cubic terms (where there are sufficient values of the 
string tension), we can measure the spread of those results to obtain a handle on 
the systematic error connected with the extrapolation. These rather large sys- 
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tematic errors could be improved if data for the string tension were available on 
finer lattices. We have systematics of 8%, 4%, 12% for SU(2), SU(3) and SU(4) 
respectively. 

2. Interpolation of the string tension. The error connected with the interpolation of 
the string tension can be easily evaluated by interpolating at values for which the 
string tension has been directly measured. Typically this error is also well below 
1%. 

3. Determination of the A parameter using the two-loop beta function. The A pa- 
rameter is determined by assuming the validity of the two-loop beta function at 
the most perturbative point. In principle we don't know at which energy scale 
the theory is well-described by two-loop perturbation theory. If the onset of the 
two-loop regime happens at lower energy scales, we could have used the second 
most perturbative point to determine Asf- This procedure gives a systematic 
error of about 4% for SU(3) and less than a percent for SU(4). On the other 
hand, it might be possible that two-loop perturbation theory fails even for our 
most perturbative point. To check that, we have used an approximation to the 
three-loop expression 

W/3i , 1 -&'(^) 



where the scheme-dependent three-loop coefficient /3f ^ has been determined in 
combination with Asf using the two most perturbative points. We found that the 
corresponding systematic error is of order 5%. Starting at slightly lower coupling 
will reduce this error significantly. 

4. Scale uncertainty due to accumulation of iterative errors This has been explained 
in Sect. El and contributes at under 1%, so can safely be neglected. 

5. Large-N extrapolation. The large-iV extrapolation uses diagrammatic predictions 
truncated to leading correction in 1/A^. The truncation error and the onset of the 
large- regime can be determined respectively by adding higher order corrections 
and by excluding points at small A^. Due to the fact that we have data only for 
three A^ values, we do not have an estimate for these errors in our calculation. 

To sum up, an estimate of the total systematic error is 15%, much greater than the 
statistical error, and almost all which comes from point 1 above. 



6 Conclusions 

In this work, we have formulated the Schrodinger functional for SU(4) lattice gauge 
theory and studied it numerically. The resulting running of the coupling seems to be 
correctly described by two-loop perturbation theory down to energy scales of the order 
of ^/a. This could be an indication of an underlying exact (5 function [29], at least 
at large A^ [30]. Our results were also used to determine A;^. Using our calculation 
and the calculations for SU(2) and SU(3) given respectively in Refs. [3] and [4], we 
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have performed a comparison with the extrapolation to ^ oo of [22], finding good 
agreement. Although the limited number of points in our study (A^=2,3,4) does not 
allow us to perform a large- extrapolation, it is reassuring that the two studies give 
compatible results, as the dominant sources of systematic errors are different in the 
two cases. In order to perform a controlled extrapolation to = oo, we are currently 
extending the calculation to = 6, 8. This will also require reducing the current 
systematic error on the value of A-p^, which mostly comes from the extrapolation to 
the continuum limit of Lmax- The smaller systematic error in SU(3) with respect to 
SU(2) and SU(4) suggests that this can be achieved by using an improved action that 
suppresses the linear term in a^/a. Finally, it would be interesting to compare the 
running of the coupling in the SF scheme to the running in the recently introduced 
interface tension scheme [31,32]. 
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